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Abstract — This study considers an orthogonal frequency- 
division multiple-access (OFDMA) -based multi-user two-way re- 
lay network where multiple mobile stations (MSs) communicate 
with a common base station (BS) via multiple relay stations (RSs). 
We study the joint optimization problem of subcarrier-pairing 
based relay-power allocation, relay selection, and subcarrier 
assignment. The problem is formulated as a mixed integer 
programming problem. By using the dual method, we propose 
an efficient algorithm to solve the problem in an asymptotically 
optimal manner. Simulation results show that the proposed 
method can improve system performance significantly over the 
conventional methods. 

Index Terms — Two-way relaying, subcarrier pairing, resource 
allocation, orthogonal frequency-division multiple-access. 



I. Introduction 

An important property of orthogonal frequency-division 
multiplexing (OFDM)-based relaying is that the frequency di- 
versity can be exploited by subcarrier pairing, which matches 
the incoming and outgoing subcarriers at the relay based on 
channel dynamics and hence improves system performance. 
In multi-user environments with orthogonal frequency-division 
multiple-access (OFDMA), subcarriers should not only be 
carefully paired at the relay but also be assigned adaptively for 
different users. If with multiple relays, it further complicates 
the problem because relay selection tightly couples with 
subcarrier pairing and assignment. Thus, subcarrier-pairing 
based resource allocation in multi-user multi-relay OFDMA 
networks is highly challenging. 

Subcarrier-pairing based resource allocation has been orig- 
inally investigated for single-user single-relay one-way com- 
munications (e.g., d, El). In particular, it is proved in 
that the ordered pairing is optimal for amplify-and-forward 
(AF) protocol. Authors in investigated separated power 
allocation and subcarrier pairing in two-way communication 
using single relay, where the power allocation is first employed 
by water-filling and then subcarriers are paired at the relay 
by a heuristic method. In 0, the subcarrier-pairing based 
joint optimization of power allocation, relay selection and 
subcarrier assignment for single-user multi-relay systems was 
studied. The subcarrier-pairing based joint optimization of 
power allocation and subcarrier-user assignment for multi-user 
single-relay scenario was studied in 0. In 0, the authors 
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Fig. 1. System model. 

studied relay-assisted bidirectional OFDMA cellular networks, 
wherein the subcarrier-pairing based joint optimization of 
bidirectional transmission mode selection, relay selection, and 
subcarrier assignment was investigated by a graph approach. 
Authors in investigated the jointly optimal channel and 
relay assignment for multi-user multi-relay two-way relay 
networks. These works 0, 0, however, did not consider 
power allocation. 

In this work, we consider an OFDMA two-way relay 
network with a common base station (BS), multiple mobile 
stations (MSs) and multiple relay stations (RSs). The downlink 
and uplink traffic for each MS is multiplexed through analog 
network coding at the RSs. We formulate a joint optimization 
problem of subcarrier-pairing based relay-power allocation, 
relay selection, and subcarrier assignment. The problem is a 
mixed integer programming problem and we solve it efficiently 
in dual domain with polynomial complexity. 

II. System Model and Problem formulation 

We consider a single-cell OFDMA two-way relay network, 
as shown in Fig. [T] with one BS, multiple MSs and RSs. 
All the MSs are assumed to be cell-edge users so that 
both the downlink and uplink traffic of each user needs to 
be relayed through one or more RSs. This assumption is 
commonly used for cellular relay networks in the literature 
(e.g., 0, 0, 0). Each RS operates in a half-duplex mode 
and relays the bi-directional traffic using AF protocol, known 
as analog network coding. In specific, the AF two-way relay 
protocol takes place in two phases ifTOl . In the first phase, also 
known as multiple-access (MAC) phase, all the MSs and the 
BS concurrently transmit signals while al the RSs listen. In 
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the second phase, known as broadcast (BC) phase, the RSs 
amplify the received signals and then forward them to the 
MSs and the BS. To avoid multi-user interference, each MS 
and RS operate in non-overlapping subcarriers in the first and 
second phases, respectively. The downlink-uplink interference 
within each user is eliminated by self-interference cancelation. 
Furthermore, the channel is assumed to be slowly time- varying 
and all the channel state information can be perfectly estimated 
and known at the BS for centralized processing. 

Let J\T = {1, 2, • • •, N} denote the set of subcarriers, JC = 
{1, 2, • • ., K} denote the set of RSs, and U = {1, 2, • • M} 
denote the set of MSs. The channel coefficients from BS b and 
MS u to RS k on subcarrier i in the MAC phase are denoted 
as hb^,i and f u ,k,i, respectively, \/u G U, k G /C, i G A/\ In 
the BC phase, the channel coefficients from RS k to BS b 
and MS u on subcarrier j are denoted as hk,b,j and fk,u,j, 
respectively, \/u G U,k G /C, j G N. Here channel reciprocity 
is used, which is valid in TDD (time-division duplex) mode. 
Along with the paths, we further denote Pb,k,% and p u ,k,i as 
the transmitted power of BS b and MS u respectively, and 
Pk,u,j as the transmitted power of RS k. Then, the sum-rate 
of uplink-downlink transmission of MS u over subcarrier pair 
(i, j) with the assistance of RS k can be expressed as 0, fTOl 

o 1 i ft , Pu,k,i\fu,k,i\ 2 Pk,u,j\hk,b,j\ 2 \ 

tt u ,k,i,j =~ log 2 1 H ^— 

2 V Pk,u,j\hk,b,jr + mu,k,i J ^ 

i 1 ! A , P6,fc,i|^b,fc,i| 2 Pfc,n,j|/fc,n,j| 2 ^ 
+ 9 io S2 I 1 H T7 |2~7 ) 5 

^ \ Pk,u,j \Jk,u,j | ~r fT^u,k,i J 

in which m w>fc ,i = 1 + P6,/c,^6,/c,i| 2 + Pn,fc,i|/n,/c,i| 2 . It can 
be proved that the sum-rate R Uyk ,ij is concave in the relay 
power Pk,u,j- 

We then introduce a set of binary variables p u ,k,i,j £ {0,1} 
for all ix, fc, z, j. When p u ,k,i,j = 1> it means that subcarrier 
z in the MAC phase is paired with subcarrier j in the BC 
phase and they are used by RS k to relay the signals of MS 
u. Otherwise, we have p u ,k,i,j = 0. These binary variables 
must satisfy the following constraints, due to the exclusive 
subcarrier assignment, 

M K N 

^2^2^2 Pu,k,i,j < 1, Vi, (2) 

u=l k=l 3 = 1 
M K N 

u=l k=l i—1 

For simplicity, we study relay-power allocation and let the 
transmit power of the BS and MSs be fixed. Each RS is subject 
to its own peak power constraint. This can be expressed as: 

M N 

^2^Pk,u,j <Pk, Vfc, (4) 

u=l j=l 

where is the peak power constraint of RS k. 

Our objective is to maximize the system total weighted 
throughput by jointly optimizing the assignment variables 
{pu,k,i,j} and the relay power variables {pk,u,j}> Mathemati- 



cally, this can be formulated as: 

M K N N 

max w u / J / ]/ ] Pu,k,i,jR u ,k,i,j(Pk,u,j) (5) 

11 u=l k=l i=l j=l 

s.t. da),®,©, 

where w u is the weight that represents the priority of MS u, 

p G Rf WxJV and p G {0, 1}^xKxATxAt are matrices with 

entries Pk,u,j and p Uj k,i,j, respectively. 

III. Dual Based Algorithm 

We first define T as the set of all possible p satisfying © 
and ©, V as the set of all possible power allocations p for the 
given p that satisfy p^ uJ > for p u ,k,i,j = 1 and p kjUij = 
for p u ,k,i,j — 0- Denote A = (Ai, A2, A^) ^ as the 
dual variables associated with the peak power constraints of 
the RSs. Then the dual function of the problem in © can be 
defined as 

g(X) = max L(p, p, A), (6) 

pev(p) 
per 

where the Lagrangian is 

M K N N 

L{p,p,\) = ^2 w u^2^2^2R u ,k,ij(Pk,uj) 

u—l k—1 i—1 j—1 

K I M N \ 

k=l \ n=lj=l J 

Computing the dual function g(A) requires us to determine 
the optimal (p, p) at the given dual vector A. In the following 
we present the derivations in detail. 

A. Optimizing the Primal Variables (p, p) for Given A 

We first find the optimal power variables p by fixing the 
binary assignment variables p. Then we search the optimal p 
by eliminating p in the objective function. Such a method has 
been commonly used in the literature (e.g., O, H, (5), fTTIl ). 

Let us rewrite L(p, p, A) as 

K U N K 
L{P,P,X) = ^2^2^2 L k,u,j(Pk,u,j) + ^ AfcPfc, (8) 
k—1 u—l j=l k—1 

where 

N 

L u ,k,j(Pu,k,j) = Wu^^Ru,k,i,j(Pk,u,j) ~ ^kPk,u,j- (9) 
i=l 

Suppose p u ,k,ij = 1 for a certain (u,k,i,j). It is easy to 
verify that L u ^j(p u ,k,j) is concave in Pk,u,j and thus the 
optimal p£ u j(Afc) can be obtained by applying the Karush- 
Kuhn-Tucker (KKT) conditions. More specifically, P%, u ,j(^k) 
is the non-negative real root of the following quartic function 

a Pk,uj + h Pk,u,j + c Pk,uj + d Pk,u,j + e = 0, (10) 

where a, 6, c, e are coefficients determined by the dual 
variables, MSs' weights, and channel gains as defined at the 
top of the next page. 
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Substituting the optimal power allocations p*(A) into © 
to eliminate the power variables, the dual function can be 
rewritten as 

M K N N K 

g(X) = max ^ ^ ^ ^ p U) jfc,i ) j^,jb ) i,j + ^ A fc P fc , 

u=l fc=l i—l j=l k—1 

where 

X u ,k,i,j = w uRu,k,i,j(Pk,uj(^k)) — ^kPk,uj(^k)' (H) 

Now we are ready to find the optimal p. In the following, we 
show that X Uykyi j defined in (fTTT) plays an important role in 
user and relay selection for occupying a subcarrier pair 

Noting the constraints © and ©, we conclude that there is 
at most one non-zero element for a given subcarrier pair 
This suggests that at most one MS and one RS can occupy 
the subcarrier pair (i, j). Based on the observation, we define 

Xij = max X Ui k,i,j, (12) 

(u*,k*)ij = arg max X ukij . (13) 
keic y ueu 

Then the dual function can be further written as 

N N K 

g(X) = max ^ ^ + ^A fc P fe . (14) 

^ 2=1 j = l fe=l 

From (fl4l) it can be seen that if subcarrier i in the MAC phase 
is paired with subcarrier j in the BC phase, then the pair should 
be used by MS u* with the help of RS fc*, i.e., the MS and 
RS with the maximum X Ujkj ij as defined in (ITTb . This can be 
interpreted from an economic perspective. Suppose each dual 
variable X k represents the power price of RS k. Then X Uj k,ij 
can be regarded as the profit of letting MS u transmitting over 
the subcarrier pair (i, j) with the help of RS k, where the 
profit is defined as the throughput revenue w u R u k % j minus 
the power cost X k p* ku -. Clearly, to maximize the system total 
profit, each potential subcarrier pair (i, j) should be assigned 
to the MS and RS that can generate the maximum sub-profit. 

The remaining problem is then to identify the optimal sub- 
carrier pairings p u *,k*,i,j- This is a standard two-dimensional 
assignment problem. The classical Hungarian method can be 
applied to find the optimal p* (A) in polynomial time. 



B. Optimizing the Dual Vector A 

After computing g(\), we now solve the standard dual 
optimization problem which is 

mm g(X) (15) 
s.t. XhO. 

Since a dual function is always convex, subgradient-based 
methods can be used to minimize g(X) with global conver- 
gence with the fact that 

M N 

AX k = P k -J2J2Ptuj^k) (16) 

U=l j=l 

is the subgradient at Afc, Vfc. In specific, denote AA^ = 
(AX±\ AAg \ AXj)), then we can update the dual vari- 
ables as A ( ' +1) = A ( "^ W } AA (0 . Here, uo^ is the dimin- 
ishing step size at the Ith iteration to guarantee the convergence 
of the subgradient method. 

C. Refinement of Power Allocation 

Having the dual point A*, we now need to determine the 
optimal solution to the primal problem Q. Due to the non-zero 
duality gap, the optimal p*(A*) and p*(A*) may not satisfy 
all the constraints ©, ©, and © in the original problem. 
To overcome this problem, we first determine the optimal 
assignment p* in dual domain, and then make a refinement 
of the power allocation to meet the power constraints in the 
primal problem. More specifically, denote A u ,k as the set of 
active subcarrier pairs assigned to MS u and RS k obtained 
from the dual problem. The problem can be written as: 

M K 

max ^2 w u ^2 Y R u,k,i,j(Pk,u,j) (17) 
u =i k=i (i,j)eA u , k 

M 

s - L p*™ - Pk > yk - (18) 

u=1 (hj)eA u , k 

Clearly, this is a convex problem. By applying KKT condi- 
tions, we can obtain the optimal p% u j which has the same 
expression as that in the dual domain. 

Finally we summarize the overall procedure of the proposed 
dual-based solution in Algorithm 1 . This algorithm is asymp- 
totically optimal when N is sufficiently large lfl2ll . 



IEEE WIRELESS COMMUNICATIONS LETTERS 



4 



Algorithm 1 Proposed algorithm for problem ([5]) 

1: initialize as a random non-negative vector, I = 0. 
2: repeat 

3: Compute X u ,k,i,j using dTTb for all (u,k,i,j) with p%, u ,j 

being the non-negative real root of ([Tol) . 
4: Obtain Xij and (u* , k*) using (fl2t and (\3l respectively for 

all then obtain optimal p*(A^) by solving ([HI) . 

5: Update A w using the subgradients AA W in dT6j); Let / «- 

1 + 1. 
6: until A converges. 

7: Set the final p as p* obtained in the dual domain and refine the 
power parameter p* by solving dTTb at the given p* . 



D. Discussion on Complexity and Proportional Fairness 

The complexity of updating the dual variables A is 0(K q ) 
(e.g., if the ellipsoid method is used, q = 2). The complexity 
in (El) and the Hungarian method are 0(MK) and (D(N 3 ), 
respectively. Combining all, the total complexity of the pro- 
posed method is 0((MK + N 3 )K q ), which is polynomial. 

If consider long-term fairness among the MSs, the weight 
of MS u at time t can be updated by w$ = 1/Tu \ e U, 
where as the accumulated rate of MS u at time t. Note 
that we can let w u = 1 for every MS for pure throughput 
maximization. 

IV. Simulation Results 

We consider a cell with 2 km radius. The RSs are uniformly 
located on a circle centered at the BS and with radius of 1 
km. The MSs are randomly but uniformly distributed in the 
outer circle as in Fig. 1. The path loss exponent is 4 and the 
standard deviation of log-normal shadowing is 5.8 dB. The 
small-scale fading is modeled by multi-path Rayleigh fading 
process. A total of 3000 independent channel realizations were 
used. Each channel realization is associated with a different 
set of node locations. We set M = 4, K = 3, and N = 32. All 
MS and the BS have the same maximum power constraints, 
so do all RSs. We set the BS and MS power to be 10 dB 
per-node and uniformly distributed among all subcarriers. 

As the benchmarks, the performance of Equal Power As- 
signment (EPA) based resource allocation and Random Re- 
source Allocation (RRA) schemes are also presented. Specif- 
ically, EPA lets p be uniformly distributed among all the 
subcarriers on each relay station and finds optimal p* as 
in Section III-A proposed algorithm. In RRA, the power is 
uniformly distributed and the subcarrier pairs and relays are 
randomly assigned. The complexity of the EPA and RRA 
schemes are 0(MK + TV 3 ) and O (TV), respectively, which 
are lower than that of the proposed algorithm. 

Fig. [2] compares the average sum-rate achieved by different 
schemes. We first observe that the proposed dual-based algo- 
rithm approaches the upper bound (the optimal dual) closely. 
This verifies the effectiveness of the dual method at large 
number of subcarriers. One also observes that the proposed 
algorithm outperforms the two benchmarks by a significant 
margin. In particular, the proposed algorithm obtains more 
than 30% and 200% throughput improvements over the EPA 
and RRA schemes, respectively. This tremendous improve- 
ment demonstrates the superiority of our proposed algorithm. 
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Fig. 2. Average sum-rate versus RS-power per-node. 

V. Conclusion 

In this work, we have studied the subcarrier-pairing based 
resource allocation in OFDMA-based two-way relay networks. 
By using the dual method, an efficient algorithm for joint 
optimization of subcarrier-pairing based relay-power alloca- 
tion, relay selection, and subcarrier assignment was proposed. 
Simulation results show that the proposed algorithm can 
significantly improve the system performance compared with 
the conventional schemes. 
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